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Abstract 


Previous  developments  of  manpower  planning  models 
involving  uses  of  goal  programming  with  embedded  Markoff 
processes  are  here  extended  in  order  to  (a)  explicitly 
comprehend  truncational  effects--e.g. ,  those  due  to 
retirement--and  (b)  allow  for  interperiod  Markoff 
transition  matrices  which  change  over  time. 


I.  Introduction 


This  paper  reports  on  some  further  developments  in  a  research  program 
being  undertaken  in  cooperation  with  OCMM  (Office  of  Civilian  Manpower 
Management).  This  part  of  the  OCMM  effort  looks  toward  improving  the  processes 
of  manpower  planning  for  the  U.  S.  Navy  by  means  of  computer  assisted  mathe¬ 
matical  models.  See,  e.g.,  [4 ] - [8 ]  and  [11]  in  the  bibliography,  which 
reference  earlier  papers  that  extended  the  "state  of  the  art"  for  modelling 
manpower  planning  by  (a)  combining  the  ideas  of  "goal  programming"  and 
Markoff  transition  processes  and  (b)  utilizing  multiple  objectives  along 
with  (c)  other  constraints  --  such  as  budgetary  limitations  --  so  that  these 
related  considerations  could  also  be  simultaneously  treated  along  with  the 
personnel  programs  (and  proposals)  to  which  they  might  be  related.—^ 

In  this  report,  these  previously  developed  models  are  further  extended 
and  refined  to  explicitly  comprehend  truncational  effects  (e.g.  those 
associated  with  retirement)  as  well  as  the  possibility  of  interperiod  Markoff 
transition  matrices  which  change  over  time.  It  is  then  also  possible,  as  we 
shall  soon  see,  to  expend  the  transformational  formulas  that  were  found  ad¬ 
vantageous  for  simplifying  and  improving  computations  for  the  original  model. 


—  See  [3  ]  for  further  detailed  development  on  the  ideas  of  goal  programming, 
multiple  objectives,  etc.  Possibilities  for  policy  evaluations  and  sensi¬ 
tivity  analysis,  as  well  as  various  byproducts  of  value  for  personnel  plan¬ 
ning  are  also  examined  in  [4],  [5],  [6],  and  [7], 


2. 


Furthermore,  although  only  one  particular  device  for  handling  truncations  is 
presented  here,  it  will  be  seen  that  others  can  also  be  used  in  addition  to 
the  "locked-in"  truncations  that  we  shall  depict  in  this  report.  (It  is  im¬ 
portant  to  observe  this  kind  of  possibility  especially  in  connect  ^n  with 
the  many  byproduct  uses  that  the  extended  model  is  designed  to  yield.) 

All  of  the  above  research  is  specifically  designed  to  tie  in  with  the 

automated  civilian  manpower  information  systems  that  the  U.  S.  Navy  is 

also 

developing.  These  models  and  their  related  developments  are,  as  will  be  seen, 
applicable  in  other  contexts  as  well. 

II.  Recapitulation  of  OCMM's  Goal  Programming  Model  for  Manpower  Planning: 

For  ease  of  reference  we  recapitulate  and  summarize  relevant  details 
•from  the  already  developed  model  which  we  have  previously  reported.  See, 
e.g.,  [^]  through  [7].  Thus  we  write 

min  2  2  ^  <E+(t)  +  E’(t)) 
k  t 

subject  to 

(1)  E  E  E  (M^^Cr)  -  E*(t)  +  E’(t)  -  gk(t) 

t-1  jsJQ 

where 

(l.l)  g.(t)  .  f  (O  -  l  (mV  ». 

and  fk(t)  is  the  ceiling  prescribed  for  the  k^  manpower  category.  The 
subtraction  from  fk(t)  in  (1.1)  yields  a  gk(t)  value  that  represents  the  net 
requirement  for  the  k^  manpower  category  after  allowance  for  the  available 


3. 


persons  remaining  in  this  category  through  the  growth  and  attrition  which 
occurred  during  the  t  previous  periods,  where 

a  =  vector  of  initial  inventory  of  personnel  in  all  job  types. 


(2) 


M  =  Markoff  transition  matrix  depicting  the  transfer  probabilities 
between  job  types 


while  (MC)  =  the  i^  row  of  MC  and  M*  =  MM. .  .M  (t  times).  Thus  J"  £  (MC)  a"j 

1  Lie  I,.  J 


k 


th 


represents  the  carry-forward  of  the  initial  states  as  they  pertain  to  k  man¬ 
power  category  comprised  of  the  collection  1^  of  job  types  i. 

The  above  model  does  not  contain  budgetary  or  other  constraints  and 
hence  we  can  look  at  it  for  purely  manpower  planning  aspects  in  a  model  where 
the  objective  is  to  minimize  a  weighted  sum  of  the  deviations  for  each  man¬ 
power  category  via 

E^(t),  *  p^gitive  or  negative  deviation,  respectively,  for 

k  manpower  category. 


(3) 


kt 


the  weight  which  is  applicable  to  k**1  manpower 
category  in  period  t. 


Using  B(t)  to  represent  an  applicable  budgetary  ceiling  for  expendi- 
T 

tures  in  period  t  and  c  (t)  ■  (c^  (t) , . . .  ,c^(t) , . . .  ,c  (t))  as  a  vector  of 
salary  rates,  we  can  extend  the  model  in  (1)  by  adjoining  constraining  budgetary 
relations  as  follows: 

min  I  E  |i. 


k  t 


kt 


[<<« +  *;<«  j 


subject  to 


t 

Z 


Z  (M*"7)  XJ(T)  -  E+(t)  +  E'(t)  -  g  (t) 
ielk  t-1  J«JQ  1  *  *  * 


DEL  c  (t)  (Mt”T) .  xJ(t)  <  B(t)  -  cT(t)  M*a 
i  t-1  jejQ  1  1 


(4) 


4. 


where  J  refers  to  the  collection  of  external  manpower  sources  and 

j  in 

xJ(t)  >  0  refers  to  the  vector  of  amounts  secured  from  source  jej  period  t. 

j  % 

As  was  noted  in  [4],  the  basic  variables  xJ(t)  only  enter  in  certain 
combinations.  Thus  it  is  possible  to  simplify  the  model  in  (4)  by  making 
certain  transformations  for  which  the  following  symbols  (and  related  definitions) 


may  be  utilized: 

(5) 


S<7)  =  E  Xj(T) 
j€Jo 
*  t-T 

Tl(t)  -  E  MC  T  §(T) 
T*1 


The  pertinent  parts  of  (4)  can  then  become  simplified  via  the  expressions 
in  (5)  --viz., 

t  t 

(6)  E  E  E  <tr'T).  xJ(t)  -  E  I  (MC‘T)  |(t)  -  E  (T|( t) )  . 

iclk  t—  1  jeJo  idk  T-l  iclk 

Furthermore,  we  can  utilize  the  definition  M°  =  I  (the  identity  matrix) 
to  obtain 

I  l  (MW1"r).  5(t)  - 

1,1^  T-l 

t  . 

(7)  -EE  M  5(t)  +  M°  §(t+l)  - 

i«Ik  T-l 

*  t-T 

-  E  M  E  M  §(t)  +  §(t+l). 

.  *Ik  T-l 


Hence,  we  also  have 

(8)  Tl(t+1)  -  MTl(c)  +  §<t+l). 

Since  the  5(t+l)  are  vectors  of  decision  variables  and  are  non-negative,  we 
can  proceed  as  in  [4]  to  replace  them  by  the  choice  vectors  T](t)  with  the 


5. 


requirement 


7](t+l)  -  M  T\<t)  >  0 


(9)  and 


T|(t)  >  0  . 


Preparatory  to  the  interpretations  and  elaborations  undertaken  in  the 

next  section,  we  can  utilize  these  developments  to  replace  the  original  model 

in  (4)  by  the  following  transformed  and  reduced  version  --  viz. , 

min  I  I  nkt  [E+(t)  +  E"(t)  ] 
k  t 


subject  to: 


(10) 


Z  \'t)  -  E^(t)  +  E‘(t)  -  gk(t) 

\ 

I  cjL(t)  ^(t)  <  B(t)  -  cT(t)  Mta 

-  (M)t  Tl(t)  +  Ts^t-vl)  >  0 
\(t),  E+<t),  E’(t)  >  0. 


III.  Interpretations  and  Further  Development: 

To  interpret  what  has  been  achieved  above,  note  first  that  T|^(t)  is 
the  ith  component  of  T|(t)  where 

(11)  T)(t)  ■  Accumulated  Inventory  from  hiring  to  period  t. 

The  vector  Tt(t)  can  be  related  to  the  recruitment  that  will  be  programmed  by 
means  of  the  vector  |(t)  and  other  relations  as  follows.  First,  the  amount 
to  be  recruited  in  period  1  is  given  by 


T|(l)  «  ?(l> 


6. 


where  5(1)  is  the  amount  to  be  recruited  for  period  1.  See  (5).  For  period  2, 
however. 


11(2)  *  M§(1)  +  1(2) 

where  5(2)  is  the  amount  to  be  recruited  in  period  2  and  M  is  the  Markoff 
matrix  which  has  as  its  elements  the  probabilities  of  movement  into  and  out 
of  the  various  jobs.  See  (2).  Oi  course,  (M)i  is  the  ith  row  of  M  and 

V 


(M)^  * 


relates  this  row  to  the  elements  of  M.  The  expression  which  relates  the 
inventory  of  hirings  in  period  3  to  the  hiring  activities  up  through  this 
period  is 

71(3)  -  M2  5d)  +  M  5(2)  +  5(3), 


and  so  on. 

This  leads  to  a  "triangular  system"  which  may  be  developed  for  convenient 
representation  and  solution  by  means  of  the  following  definitions 


I  o' 

I  0 

'V 

H 

ft _ 

■ 

l. 

M  1 

m  a 

X+M  Ij 

(12)  where 

1  •  Identity  matrix 
X  •  -M. 


Thus,  for  the  2  period  case: 

I  0 


I 

M 


and 


*71(1)1 

.W>J 


7. 


’5(1)' 

gg 

I 

on 

-1 

H 

pr 

L 

M 

_TK2>J 

1 

« 

0 

’11(1)' 

a 

M 

I. 

.H(2). 

since  —  see  (1?)  —  the  indicated  inverse  exists. 
Similarly  for  3  periods 


5d)' 

o 

o 

-1 

11(1)' 

"l  0  o' 

il(D' 

5(2) 

m 

M  I  0 

11(2) 

S 

-M  1  0 

11(2) 

5(3)_ 

M2  M  I 

11(3) 

- v 

M 

a: 

1 

o 

_ * _ 1 

11(3) 

?(D  -  “fid) 

5(2)  -  11(2)  -  MT1(1) 

50)  -  11(3)  -  MTl(2), 

and  so  on. 

More  generally,  we  can  write 

(13)  5(0  ■  H(0  *  Ml}(t*l). 

Pausing  and  reflecting  on  this  structure  we  see  that  it  can  be  significantly 
generalised  to  the  situation  where  the  Markoff  matrices  say  differ  from  period 

to  period.  Thus  if  li/t)  is  accumulating  through  period  t  we  can  write 

(14)  5'0  •  m(0  *  l!(t-i) 

I 

where  t  j  is  the  probability  matrix  applying  to  the  transitions  from  period 
t*l  to  t.  'Note  that  the  inversion  is  still  strai£htforvard  in  any  case  and 
requires  no  division  operations,  whether  M  is  non-sing  :lar  or  not.) 


8. 


IV.  Adjustments  for  Variations  in  Proportions  Eligible  for  Retirement: 

Part  One:  Estimation  and  Adjustment  Formulas 

We  now  develop  a  procedure  for  further  extending  the  OCMM  goal  pro¬ 
gramming  model--or  any  other  such  model--and  its  application  by  developing 
certain  formulas  which  make  better  allowance  for  potential  retirements.  The 
nature  of  the  approach  to  be  igges*-ed  will  be  sufficiently  clear,  we  think, 
if  we  proceed  by  reference  to  a  five-period  example  by  means  of  which  we  can 
provide  more  specific  and  concrete  interpretations. 

Those  who  are,  or  will  be,  eligible  to  retire  in  the  next  five  years 
are  represented  by  the  vector  b, 

b  -  b5  +  b4  +  b3  +  b2  +  b1  +  b° 

where 

b^  *  those  first  eligible  to  retire  in  period  5 
4 

b  «  those  first  eligible  to  retire  in  period  4 
3 

(15)  b  *  those  first  eligible  to  retire  in  period  3 

2 

b  ■  those  first  eligible  to  retire  in  period  2 
b*  *  those  first  eligible  to  retire  in  period  1 
b^  *  those  eligible  to  retire  in  period  0,  but  did  not  do  so. 
Next,  let  a  *  a  +  b  where  a  *  total  "on  board"  at  t  =  0,  a  =  those  on  board 
not  eligible  tj  retire  in  the  next  five  years  and  b  *  those  on  board  who  are 
eligible  to  retire  in  the  next  five  years. 

Now  we  introduce  the  following  expressions  and  definitions: 

«  «  C0/B0  ,  i  -  0  -  C0/B0> 


where 


9. 


(16) 


Cq  *  Personnel  remaining  from  those  eligible  to  retire  in  some 
base  period 

B  *  Total  eligible  to  retire  in  some  base  period. 


For  those  who  are  eligible  to  retire  in  the  next  five  periods,  we  shall 

assume  that  either  they  will  retire  or  that  they  will  remain  in  the  same  job. 
our  formulas 

Then  we  develop  for  what  we  shall  call  a  five-period  "stretch".  This  is  done 
A 

by  period,  as  follows,  where  the  formula  for  estimating  those  who  will  not 
retire  is  shown  first  in  each  of  the  five  periods: 

Period  1: 

a  (b1  4  b°) 
a  (b1  +  b°) 

Period  2: 

-  2  -2  1 

+  Qf  (b  +  b  ) 

orb2  4  a<5(b*  +  b°) 

(17)  Period  3: 

ab3  +  a  [ab2  4  a^b1  4  b°)] 

ab3  +  a  [ab2  4-  a2(b1  4-  b°)  ] 

Period  4: 

ab4  4-  a  [ab3  4-  a2b2  4-  i3(b1  4  b°)  ] 

ab4  4-  a  [ob3  4  a2b2  4  c?3(b1  4  b0)] 

Period  5: 

ab3  4  a  [ab  ^  4  a2b3  4  a3b2  4  a4(b^  4  b°)  ]. 

ob5  4  a  [ab4  4  a2b3  4  a3b2  4  a4  (b1  4  b°)  ] 

Finally,  it  also  seems  fair  to  assume  that  those  who  are  recruited  during 

this  five-period  stretch  will  not  be  eligible  to  retire  during  this  stretch. 


in 


UIL  T wo:  Goal  Adjustments 

•  d*  extended. model  nay.be  secured  from. (10)  by  using  the  precedi t\ 


and  a* sumptions  to  adjust  for  the  indicated  retirements.  These 
"retirement  adjustment’1  corrections  will  be  reflected,  of  course,  in  the 
budget  constraints  as  veil  as  in  the  goals  and  objectives. 

da  was  observed  in  section  III,  .the  original  model  utilized  the  matrix, 
df,  hs  prajaevt  ell  of  the  initial  taaepoax..  plus  maur owe .•  '.lire  '  di  v •  t' 
~wotr»tch".-^  hfc  now  propose  to  utilize  our  preceding  development,  however,  u. 
encompass  the  situations  involving  retirement  projections. 

A  k 

Suit,  if  w*  use  Ti(fc)  to  represent  the  cumulative  manpower  vector  through 
period  t,  then  we  can  represent  as  a  sum  of  two  vectors  —  viz., 

(18)  1}(t)  -  T,l(t)  +  H2(t) 


where  T]  (t)  represents  the  manpower  cumulation  not  affected  by  retirement  and 
2 

T)  (i)  represents  the  cumulation  affected  by  retirement.  Furthermore,  by  an 
evident  extension  of  (13)  joined  to  the  development  in  Part  One  of  this  sec  tier 


(19) 


TlX(t)  -  MTlV-l)  +  §(t) 
*n2(e)  -  aTl2(t-l)  +  ab*. 


Specifically, 


(19.1) 


Tj2(l)  «  a  (b1  +  b°) 

Tl2(2)  «  ob2  +  aTf (1) 

-.2  .  -2,,  1  .  .0. 

■  orb  +  o?  (b  +  b  ) 


as  required  for  period  (2)  --  see  (17)  — 
Now  we  let 

(20)  a^  *  those  in  job  i  initially 


1/ 


E.g.,  the  five-period  stretch  which  we 


and  so  on. 

(at  t  *=  0) 
are  using. 


11. 


and 


l 

I 

f 

i 

i 


(21)  *  at  +  b^, 

so  that,  as  in  the  development  which  immediately  follows  (15),  we  have  a^  repre¬ 
senting  those  not  eligible  for  retirement  in  the  five-period  stretch  covered  by^' 


|  (22)  b  =  b°  +  b1  +  ...  +  b5. 

f  The  vector  of  personnel  in  period  1  coming  from  within  is 

\ 

(23)  M  a  +  a  (b°  +  b1)  =  M  a  +  ^(l) , 

via  the  first  expression  in  (19.1).  Similarly,  the  vector  of  personnel  in 
period  2  coming  from  within  is 

(24)  M2  a  +  a  b2  +  a2  (b°  +  b1)  =  M2  a  +  Tl2(2), 


via  the  second  expression  in  (19.1). 

Now  the  total  on  board  in  period  1  equals 

(25)  M  a  +  T|l(l)  +  T]2(l)  +  §(1) 

where  §(1)  represents  hiring  from  outside  while  M  a  provides  the  transitions 
from  the  state  given  by  a,  the  remnant  from  t*  0.  The  vector  of  pe  sonnel 
from  within  at  period  two  is 

(26.1)  M2  a  +  M  §(1)  4  T!2(2) 


which,  in  turn,  forms  a  part  of  the  total  period-two  personnel 
(26.2)  M2a  +  M§(1)  +  T)2(2)  4  §(2). 


-^See  (15). 


12. 


By  an  evident  extension  we  can  now  write 


(27) 


Mfca  +  r  Mt-T  §(t)  +  T]2(t) 
T=1 


as  the  expression  for  the  vector  of  total  personnel  in  period  t.  By  virtue 


of  the  definition  of  T)  (t)  in  (18), 

t 


t-1 


(28)  T)X(t)  =  2  Mt'T§(T)  s  §(t)  +  2  Mt_T  |(t), 

T*1  T= 1 


But, 

(29) 


t-1  t-1  .  ..  - 

2  M  §(T)  =  M  2  M(  §(t)  =  T1  (t-1). 

T*1  T— 1 


Thus,  via  (28)  and  (29), 
(30) 


Tll(t)  =  5(t)  +  TJ1  (t-1) . 


V.  Adjustment  and  Extension  of  the  Model: 

Now  turning  to  the  model,  (10),  the  objective  function  we  wish  to 
minimize  is 

•  „ 

2  2  ^  (E+(t)  +  E*(t)) 

k  t 

subject  to  goal  (e.g.,  manpower  cei.lngs)  and  budget  constraints, 

tH 

where  Ujct  *  weights  associated  with  the  k  manpower  ceiling  in  period  t 

til 

Ek(t)  *  excess  above  k  goal  in  period  t 

Ek(t)  *  deficiency  below  k^  goal  in  period  t. 

Both  Ek(t)  and  Ek(t)  will  not  appear  in  the  solution  at  the  same  time  since 
a  given  goal  discrepancy  can  only,  be  in  one  direction  at  any  given  time 
(e.g.,  if  E^(t)  is  in  the  solution  then  E^(t)  will  equal  zero). 


13. 


Developing  the  goal  constraints,  in  which  the  a  E*(t)  and  E^(t) 
th 

variables  appear,  the  k  goal  is  defined  as 

£  [Total  personnel  in  period  t],  (t)  +  E  (t)  “  f  (t) 

ieL 


where 


fk(t>  =  ktn  manpower  goal  (or  ceiling)  in  period  t. 


Substituting  the  expression  for  total  personnel  in  period  t  we  have 


I  [<Mt5)1  +  At)  Tl (c)  ]  -  E*  (t)  +  e'  It)  -  f  (t) 
UIk  ■  '  ' 


Rearranging, the  goal  constraint  is 

I  -  st(t)  +  E‘(t)  •  fk(t)  -  r  (Hta>1  -  E  At)  -  8k(t) 

iclk  UIk 

where  g^(t)  Is  the  k^  manpower  requirement  obtained  from  the  goal,  net 

after  allowing  for  the  available  persons  remaining  in  this  category  from  the 

initial  population. 

The  budget  constraints  are  defined  as 

T 

[Cost  vector  in  period  t]  [Total  personnel  vector  in  period  t] 

<  B(t) 

where 

B(t)  ■  total  dollar  budget  for  personnel  in  period  t 


14. 


and  the  T  superscript  stands  for  "transfers",  as  in  (4).  Substituting  the 
expressions  for  the  cost  vector  in  period  t  and  for  the  total  personnel  vector 
in  period  t  we  obtain 

cT(t)[Mti  +  V(t)  +  T]2(t)  ]  <  B(t) 


where 


c  (t)  £  [c1(t), . . . ,cn(t) 1. 


Rearranging, we  obtain  the  budget  constraint 

cT(t)  Vet)  <  B(t)  -  cT(t)  M ll  -  cT(t)  T)2(t) . 

Thus,  in  conclusion,  the  model  (10)  can  now  be  replaced  by 

min  E  E  (E+(t)  +  E"(t)] 

k  t 


Subject  to: 


lei 


Z  TlJ(t)  -  E+(t)  +  E”(t)  -  gk(t) 


(30) 


cT(t)  V(0  <  B(t)  -  cT(t)  MCa  -  cT(t)  V(0 


V(t)  >  H  T|l(t-1) 


Ev,  E"  >0. 


Should  one  wish  to  go  atill  further,  e.g.,  to  also  permit  the  transition 
matrix  to  change  from  period  to  period,  all  that  would  change  in  the  above 
expressions  would  be  to  replace  MC  by  Kt-l)  ***  Kl)  *n<*  M  By 

See  (14)  ff. 
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